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ABSTRACT 

Observational and theoretical studies of extragalactic radio sources have suggested 
that an inhomogeneous environment may be responsible for observed arm length asym- 
metries of jets and the properties of extended emission line regions in high redshift 
radio galaxies. We perform 3D hydrodynamic simulations of the interaction of a pow- 
erful extragalactic bipolar jet with a disc-shaped clumpy interstellar medium of log- 
normal density distribution and analyze the asymmetry. Furthermore, we compute the 
relation between jet asymmetry and the ISM properties by means of Monte Carlo sim- 
ulations based on a ID propagation model for the jet through the dense medium. We 
find that the properties of the ISM can be related to a probability distribution of jet 
arm length asymmetries: Disc density and height are found to have the largest effect 
on the asymmetry for realistic parameter ranges, while the Fourier energy spectrum 
of the ISM and turbulent Mach number only have a smaller effect. The hydrodynamic 
simulations show that asymmetries generally may be even larger than expected from 
the ID model due to the complex interaction of the jet and its bow shock with gaseous 
clumps, which goes much beyond simple energy disposal. From our results, observed 
asymmetries of medium-sized local radio galaxies may be explained by gas masses of 
10 9 to 10 10 Mq in massive elliptical galaxies. Furthermore, the simulations provide a 
theoretical basis for the observed correlation that emission line nebulae are generally 
found to be brighter o n the side of the shorter lobe in high redshift radio galaxies 
([McCarthy et al.lll99lf ). This interaction of jets with the cold gas phase suggests that 
star formation in evolving high redshift galaxies may be affected considerably by jet 
activity. 

Key words: galaxies: jets, galaxies: ISM, methods: numerical, hydrodynamics ISM: 
structure 



1 INTRODUCTION 

Asymmetries in double radio sources have been a matter of 
discus sion since more than 40 years, when iRvle & Longairl 
(1967) examined a conjectured relativist ic propagation. This 
discussion has included asymmetries both in brightness 
and in length (e.g. of the lobes or the core-hotspot dis- 
tance). The underlying assumption mostly has been that 
jets are intrinsically symmetric objects as they are launched 
near the central supermassive black holes. While relativis- 
tic effects may cause brightness asymmetries by Doppler 
beaming ([Wardle &; Aaronl Il997h and length asy mmetries 
due to finite light t ravel times (exam i ned by IScheuei 



199 j : lArshakian fc Longairl l200d : Ijevakumar et al] 12005) 



1995 ), several studies (iPedeltv et al.lll989l: iMcCarthy et al 
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have suggested that an inhomogeneous environment has a 
stronger impac t on the obser ved asymmetries. This is in 
agreement with IScheuerl (|l995l ) finding only a slow jet head 
propagation to be compatible with observations. In contrast 
to the local universe, where low gas mass fractions are gen- 
erally found in radio galaxies, galaxies at high redshift show 
gas masses of t he same order as t he stellar mass of the 
hosting galaxy ([Tacconi et al.l I201Q.). Such large amounts 
of gas in and around galaxies may have considerable ef- 
fects on the jets. The interaction of jets with the clumpy 
interstellar medium (ISM) and gas accreted onto the galaxy 
may have a major impact on the formation of stars in the 
galaxies and the galaxies' energy budget. These questions 
go much beyond the merely morphological properties of the 
radio galaxies, but are clearly linked with the kinematics, 
dynamics and evolution of the galaxy as a whole. Obser- 
vations of high-redshift radio galaxies commonly show the 
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Figure 1. 3D density volume rendering of the jet pushing through 
the dense gaseous disc at t = 5.4 Myr, box size 32 x 32 x 16 kpc 
(cut in the plane of the jet). Low density is indicated by red, high 
density by blue. 



"alignment effect" of having extended and sometimes very 
massive emission - line nebulae aligned w ith the radio source 
(|McCarthvl[l993l : iNesvadba et alJl2QQ8h . These nebulae al- 
most always are brighter in [Oil] emission on the side of 
shorter jet. This suggests that there may be vivid interac- 
tion going on: assuming the brighter shock-excited gas traces 
the location of stronger interaction of the jet with its en- 
virons, this may be an indication of the jet being respon- 
sible for the formation of t he emission line nebulae either 
direc t ly or via its backflow (Gaibler , Krause fc Camenzindl 
l2009l : iGaibler fc Camenzindll2010h . 

Aiming towards understanding the interaction of jets 
with a complex multiphase environment at high redshift, we 
have performed 3D hydrodynamic simulations of a pair of 
jets with a surrounding clumpy gaseous disc. In this paper, 
we present first results from the simulations, concentrating 
on jet asymmetries caused by the clumpy ISM. We present 
results of Monte Carlo simulations based on analytical ap- 
proximations of the jet propagation through a clumpy ISM 
aside with the full hydrodynamical simulations to examine 
the link between the gaseous disc properties with the ex- 
pected jet length asymmetries. A more detailed examination 
of the simulations with respect to star formation and feed- 
back energetics as well as the kinematics of the dense gas 
will be presented in a subsequent paper. 



2 THE MODEL 

We set up a clumpy gaseous disc contained within a homoge- 
neous hot atmosphere (Fig.[T]). The atmosphere has density 
of p a = 0.05 m p cm -3 and a temperature of T a = 1.15 x 10 7 
K. The log-normal density field of the disc is constructed 
from a normal-distributed (mean 0, standard deviation 1) 
and clumpy field /(x), constructed in Fourier space, by 

p(x) = p exp{cr/(x)} exp(— R/Ro) sech 2 (h/ho) , (1) 
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Figure 2. Energy spectrum E(k) of the density field at the initial 
conditions. The dotted line shows a power law E(k) oc fc- 5 / 3 for 
comparison. The bump at k ~ 40 kpc -1 is a sampling effect since 
the disc cells are on a grid a factor of 2 coarser than the finest 
grid. 



with the central median disc density p — 10 m p cm -3 and 
the disc scales Ro = 5 kpc and ho = 1.5 kpc (R, h: cylindri- 
cal radius and vertical height, respectively), using a temper- 
ature of 10 4 K. The disc values are only set where p(x) > p a , 
and a strict cutoff is applied at R = 16 kpc and h = 6 kpc. 
The total mass of the gaseous disc is 1.5 x 10 n M©. For any 
small volume, the disc densities show a log-normal distribu- 
tion, but on the global disc scale, the imposed radial and ver- 
tical profiles increase the contribution of lower densities. The 
standard deviation a of the disc's log-density corresponds to 
turbu lence of Mach number 5 (~ 80 km s _1 ) (Kr itsuk et al.l 
l2007l ), a value similar to large evo lving discs at high red- 
shift dForster Schreiber et al7l2009h . The two-point struc- 
ture of the density field is described by the Fourier energy 
spectrum E{k) (Fig. [2]), which follows E(k) oc k~ 5/3 for 
large wave numbers k, but is intentionally damped towards 
smaller wave numbers, avoiding large inhomogeneities on 
scales I > ho. The initial configuration is evolved for 2 Myr, 
and after this initial relaxation phase a bipolar (back-to- 
back) jet is introduced in the centre of the disc. It is realized 
by a cylindrical orifice of jet plasma (pj = 5 x 10 -5 m p cm -3 , 
Vj = 0.8 c) with the same pressure as the hot atmosphere 
and with a radius of rj = 0.4 kpc and an initial length in 
both directions of 3 rj , respectively. The kinetic power of the 
jet is Lkin = 5.5 x 10 45 erg s _1 . 

This setup is evolved by the RAMSES 3.0 code (iTevssierl 

l2002h . a non-relativistic second-order Godunov-type adap- 
tive mesh refinement (AMR) code. The total computational 
domain extends over a cubic box of size 128 kpc, with a 
coarse grid of 1 kpc resolution and adaptive refinement down 
to 62.5 pc (effective resolution 2048 3 ). We use the HLLC 
solver and the MonCen slope limiter, refine on the entire 
disc and on 10 per cent gradients in density or pressure or 
speeds above 0.1 c, and assume 7 = 5/3 for all phases. This 
enables us to study virtually the entire jet-disc system at 
the finest resolution. We include the important effect of ra- 
diative cooli ng (as imp lemented in the code) for a metallicity 
Z = 0.5Z© terbll2008h . but exclude gravity and stabilize the 
disc by imposing a minimum temperature of 10 4 K since we 
do not include stellar feedback, making the disc sufficiently 
stable over the simulated time scale. 
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3 SIMULATION RESULTS 

The jet inflates two strongly over-pressured cavities filled 
with diffuse jet plasma in the central region of the disc and 
forms a disc-like structure of strongly compressed gas in the 
midplane due to the synchronous start of the two jets (the 
evolution is shown in Fig. [3]). The gas in the midplane of 
the disc at larger radii is largely unaffected by the jet due 
to the high densities there, while the cavities expand in ver- 
tical direction and form a pressure- driven blastwave in the 
outer layers of the disc and the surrounding halo gas. The 
jet beam, however, is still restrained by dense clumps and 
filaments. This early phase corresponds to the "flood and 
ch annel" and the "energy - driven bubble" phases described 
bv ISutherland & Bicknelll (|2007h . At this stage, a remark- 
able asymmetry between the two jets develops. The right- 
hand jet is only slightly deflected by a dense clump in its 
way, which it pushes sideways and considerably ablates, and 
then quickly pushes forward to the expanding bubble which 
it pierces at t ~ 5 Myr ("jet breakout phase"). The left-hand 
jet, however, is still restrained by dense clumps at somewhat 
larger radii (middle plots in Fig. [3j), causing a delay in its 
propagation with respect to the opposite jet. It reaches the 
expanding blastwave shell much later, eventually piercing it 
at t ~ 8 Myr. Already in the earliest (blastwave) phase, the 
jet's environment becomes asymmetric since the blastwave 



is able to clear more of the environment on the right side 
than on the left side, where a group of dense clumps makes 
this blastwave clearing less effective. It is also worth noting 
that the restraining clouds are not simply the ones located 
there already in the initial conditions, but they are generally 
moved and reshaped by advection and compression caused 
by the blastwave and the shocked beam plasma ( "second or- 
der" restraints). In case of the left jet, one of these clouds was 
just in the middle of the jet's path. During this phase, the jet 
plasma exhibits a complex velocity structure, forming back- 
flowing vortices in some regions and forward- flowing streams 
and plumes in otherfl Once the jets have pierced the ex- 
panding shell, they show the classic morphology known from 
jet simulations, propagating faster through the lower density 
environment, although the lateral expansion of the cocoons 
is still affected by the pierced shell. 

To examine, whether this asymmetry could have been 
expected from the clumpy disc set up in the initial con- 
ditions, we approximate the propagation of the jet head 
by one-d imensional mo mentum balance for very underdense 
jets (e.g. lKrausdl2QQ3h 

« (Pj/P) 1/2 (2) 
1 Movies available online make this more evident. 
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Figure 4. Momentum-balance estimate of jet propagation (thin 
lines, derived from initial conditions) and actual measurement 
(thick) for the left (solid) and right (dashed) jet. 

and find a propagation time out to a distance x of 

t=^~ r P (x') i/2 dx' . (3) 

Fig.[Hshows the analytical estimate along with the mea- 
sured position of the jet head based on the initial density 
field. The right jet is described quite well by the estimate, 
while the left one needs considerable additional time to 
break out compared to the expected behaviour (2.7 Myr vs. 
0.7 Myr expected from momentum balance). We attribute 
this to the mentioned "second-order" restraints - changes of 
the surrounding ISM caused by the jet activity itself. They 
are small at early times, but become important later once 
the jet has propagated through most of the disc and hence 
changed its environment significantly from the initial con- 
figuration. The asymmetry of the jets is also understandable 
in light of the initial disc mass in the jet beam's path, which 
is about 2 x 10 7 M for the right jet and 4 x 10 7 M for the 
left jet. 



4 MONTE CARLO SIMULATIONS OF 
DISC-INDUCED ASYMMETRIES 

We here investigate how the disc properties affect the distri- 
bution of asymmetries in a statistical sense, by generating 
a large ensemble of discs. Which disc leads to which typical 
asymmetries? Once the jets have broken out of the disc or 
any other clumpy ISM, the major effect on the jet morphol- 
ogy will be a propagation delay between the jets on both 
sides: in general, one jet will have propagated already some- 
what further than the other. Correspondingly, an absolute 
length asymmetry is expected rather than a relative length 
asymmetry as in the case of asymmet ries caused by light 
travel time ([Arshakian & Longairll2000h . Due to the compu- 
tational demands, a statistical treatment of this propagation 
delay by means of 3D hydrodynamic simulations is hardly 
possible. However, despite the considerable differences be- 
tween the ID estimate and the true propagation caused by 
individual clumps, simple momentum balance works well 
enough to obtain first order estimates of the effects of a 
clumpy interstellar medium on the jet propagation and the 
resulting arm length asymmetries, as well as the impact of 
disc parameters on it. Clearly, one has to keep in mind that 




time delay between jets [Myr] 

Figure 5. Probability density function derived from a Monte 
Carlo run based on the ID momentum balance estimate for p = 10 
m p cm -3 , a = 1, ho = 1.5 kpc and the energy spectrum used for 
the hydrodynamic simulation. The locations of a set of percentiles 
are indicated by lines attached to the upper coordinate axis. 

the "second-order" effects may in general cause even larger 
ranges of asymmetries than expected from a ID estimate. 
We note that while we have assumed a fully collimated, non- 
relativist ic jet for our estimates, the Monte- Carlo approach 
in combination with the one-dimensional momentum bal- 
ance could also be used for a conical and/or relativist ic jet, 
which may be more realistic for modelling the innermost 
regions, and any other ISM density structure. 

To examine the impact of disc parameters on the asym- 
metries caused by the disc, we performed Monte Carlo simu- 
lations of the jet propagation using the momentum balance 
estimate above. We varied the median value p and the log- 
normal distribution width a of the disc density, the disc 
scale height ho, the power spectrum shape and its slope at 
large wave numbers (Figs. [5] and [6}. The asymmetries are de- 
scribed by a time delay - the time difference between both 
jets once the disc has been left behind. A typical jet head 
propagation speed in the ambient medium would allow to 
relate this to an asymmetry in length. Since a homogeneous 
distribution is used for the diffuse ambient gas, no further 
asymmetries are caused at a later time within this statistical 
model and the asymmetries are entirely due to the clumpy 
disc. Fig.[5]shows the probability density for the fiducial val- 
ues of the hydrodynamic simulation. The time asymmetry 
of 0.7 Myr seen in Fig. 2] for the ID propagation model of 
our hydrodynamical setup lies well between the 50 and 75 
percentiles (cumulative probability) and hence can be con- 
sidered a "one-<7 event" . In the following, we concentrate on 
the effects of the mentioned parameters on the 50, 75 and 
90 percentiles of the delay distribution. Only one parameter 
is varied at a time while the others have their fiducial value. 

The mean density of the disc has a large impact on the 
jet asymmetries. For median densities much larger than the 
diffuse background gas, the delays vary oc p 1 / 2 , which is un- 
derstandable from 0; clearly, this also translates directly 
into changes of the total disc mass. The density variation 
within the disc, described by the log-normal distribution 
width parameter a, can be related to the turbulent Mach 
number M of the disc by 

a 2 = In (1 + b 2 M 2 ) (4) 

with b w 0.26 dKritsuk et al.1 l2007h . For the effects of a, 
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Figure 6. Impact of ISM disc properties on the distribution of the jet propagation asymmetry derived from the Monte Carlo simulations, 
plotted for the 50, 75 and 90 percentiles. Top left: disc density p, top right: density variations a (related to ISM turbulence), bottom left: 
disc scale height ho, bottom right: Fourier energy spectrum shape and slope a for large wave numbers s of jets in terms of time delays for 
constant (for the latter, thin lines are for pure power law energy spectra, thick lines are for the case with damped small wave numbers 
as in the simulation initial conditions). 



it is important to consider that for a log-normal distribu- 
tion, the median value p relates to the mean density as 
(p) — p exp(cr 2 /2), resulting in an increasing mean density 
and hence total disc mass M = AnRoho{p) for increasing a 
and constant p. To avoid effects of changing density, the top- 
right plot in Fig. [6] displays the changes in delay for constant 
mean density (i.e. constant disc mass; the median density 
changes along with a). It shows an increase in delay until 
a ~ 2 where it then decreases again. While the former is a 
result of a more and more inhomogenous disc density distri- 
bution, a very wide range of densities (for a highly turbulent 
medium with large a) results in a small median density and 
only small filling factors taken by gas of densities larger than 
the background medium value and thus lower probabilities 
for clumps being in the jet's path (for a constant mean den- 
sity, the volume filling factor of dense clumps in a log-normal 
density distribution decreases oc exp(— a 2 /8)/a for increas- 
ing a). For a constant median density, however, the delay 
times would increase even stronger for a > 2 due to the enor- 
mously increasing total disc masses. Interestingly, the delays 
for the wide range of turbulent Mach numbers 2 to 30 (cor- 
responding to 0.5 < a < 2) are only weakly dependent on 
the exact value of turbulence. The disc scale height ho shows 
a near-linear behaviour of the delays with some deviations 



for values ho similar to the initial jet length, which are an 
artefact of not modelling the earliest phase of jet propaga- 
tion. For zero initial jet length, the "convergence point" at 
ho ~ 0.5 kpc is located at the origin. In contrast to the pre- 
vious dependencies, the effect of the shape and large wave 
number powerlaw slope a of the Fourier energy spectrum 
E(k) is quite small since even without damping of small 
wave numbers these are suppressed by the imposed vertical 
profile scale ho. Furthermore the effects of small scale per- 
turbations (stronger for natter energy spectra, a — ► 0) are 
increasingly averaged out over the length scale of ho. We 
conclude from the Monte Carlo simulations that the distri- 
bution of disc-induced asymmetries mostly depends on the 
average disc density or total mass and to a weaker extent 
on the density variations within and scale height of the disc 
for plausible values of the latter two. 



5 DISCUSSION & CONCLUSION 

For a given propagation delay and once the jet has propa- 
gated beyond the vertical disc/ISM boundaries, the observed 
length asymmetry will in general be influenced by the den- 
sity profile of the ambient medium, which itself governs the 
propagation of the hotspots and lobes. Double radio sources 
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Fig ure 7. Arm len gth asymmetries for the chosen sample from 
iMullin et al.l |2008): large (top) and only small (bottom) sources. 
Radio galaxies are represented by open circles, quasars by filled 
circles. 



in this stage are expected to exhibit a length asymmetry 
Al = Vh At dependent on the propagation delay At due to 
the clumpy medium and the typical jet propagation speed 
Vh outside. The absolute length asymmetries due to the ISM 
should be roughly independent of the source size if the den- 
sity profile doesn't lead to a large change in vh> We have 
compared this to the observed asymmetries in radio galax- 
ies, using the 3CRR source s with with FR II morphology 
and z < 1 from lMullin et all (|2008j flfor a ACDM cosmology 
with Qa = 0.7, n m = 0.3 and h — 0.7. Of these, only sources 
with clearly identified hotspots on opposite sides are consid- 
ered (88 sources) and projection effects are not accounted 
for here (differences should be less than 50 per cent for most 
radio galaxies) . Although the relative asymmetries (with re- 
spect to the total source s ize) decrease for larger sources 
( Ars hakian fc Long air 2000), the absolute arm length asym- 
metries increase much. For the smaller sources (sum of LLS 
< 100 kpc, 20 sources), however, the scatter does not vary 
strongly with size and we argue that asymmetries caused 
by a clumpy ISM cannot be larger than the scatter seen on 
these scales (mostly between and ~ 10 kpc). The corre- 
sponding percentiles are 7.6 kpc (50 per cent), 10.2 kpc (75 
per cent) and 11.8 kpc (90 per cent). Based on a fiducial jet 
propagation speed of 0.05 c, the corresponding percentiles 
for the propagation delay are 0.50, 0.66 and 0.77 Myr. Com- 
parison of these values with our Monte Carlo runs now rises 
on the question, how typical situations as in our hydrody- 
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namical simulation with asymmetries larger than expected 
are. If this is the rather typical case, and actual asymmetries 
may be factors of a few larger than the ID estimates, the 
observed typical asymmetries could be described by asym- 
metries in the Monte Carlo runs of 0.1-0.2 Myr, suggesting 
disc densities are a factor of > 10 smaller than in our chosen 
set of parameters (with other combinations of parameters 
possible, or course). Although a more detailed comparison 
with observed data is beyond the scope of this paper, the 
corresponding masses for a dense clumpy medium of 10 9 
to 1O 1O M may be compatible with actually measured ga s 
masses (as e.g. in lCatinella et al.ll2010l ; lEmonts et al.ll2010h . 
However, if the asymmetric clearing of the central region by 
the initial blastwave and clumps remaining in the jet's path 
are rather rare, the observed asymmetries require disc densi- 
ties almost as high as in our hydrodynamical setup even for 
the local universe, with total masses of 1O 1O M0 and more. 
To us, the first case seems more probable since the asym- 
metric blastwave simply is caused by the asymmetric initial 
density and hence no exceptional behaviour. 

While the asymmetries of small radio sources (below 
100 kpc) may be explained by a clumpy ISM, this is not 
true for the larger sources with larger absolute length asym- 
metries. Although propagation of delayed jets through a 
declining density profile can result in increasing absolute 
asymmetries, it is hardly possible that they grow by a fac- 
tor of 10 or more. These large sources might actually become 
more asymmetric by a strongly asymmetric am bient density 
profile at large scales ([Jevakumar et al.l 12005) or due to an 
unstable jet propagation on one side possibly caused by in- 
stabilies or obstructing clumps in one jet. For the largest 
sources with large propagation speeds caused by a very 
low density environment, light travel time effects may actu- 
ally play a significant role even if it doesn't for the smaller 
sou rces. 

iMcCarthv et al.l {l991) have found that in almost all 
cases of their sample, that the extended optical line emission 
in [Oil] is brightest on the side of the shorter lobe, indicating 
that environmental effects cause the length asymmetries in 
powerful radio sources. Our 3D hydrodynamical simulations 
are now able to provide theoretical support for this scenario. 
Dense clouds or atomic or molecular gas in the ISM are able 
to restrain a propagating jet for some time due to their large 
inertia, causing asymmetries in the observed radio source. At 
the same time, the impact of the jet has considerable effects 
on the clouds themselves: depending on the exact location 
and cloud properties, they are accelerated, shock-ionized and 
ablated by the jet's ram pressure and compressed by the high 
pressure in the jet coco on. Cooling by line em ission stabilizes 
them additionally (cf. iMellema et al.l l2QQ2h and results in 
bright emission at the location of the strongest interaction 
with the jet. Fig. [8] shows the density of the group of clumps 
responsible for the large asymmetry in the simulation and 
the strongly increased pressure (hence strongest expected 
shock ionization) at the impact location of the jet on the 
clumps. The long-term survival of the clouds, however, does 
not only depend on their direct interaction with the jet and 
its bow shock, but also on the ir long exposure to the fas t 
and turbulent cocoon plasma (|Krause fe Alexander] l2007l \ 
which may destroy but also create emission line clouds. 

While radio sources in the local universe can grow more 
easily due to small gas fractions and deposit their energy 
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Figure 8. 3D Volume rendering of the interaction of jet with a 
group of clumps that caused the additional delay in the hydro- 
dynamic simulation compared to the ID propagation estimate. 
Top: density, with jet plasma tracer overlaid in yellow. High den- 
sities are blue, lower densities green. Bottom: pressure, with high 
pressures in red and low pressures in blue. The boxes are 8.8 kpc 
wide and 1.8 kpc thick, t = 4.84 Myr. 



in or near the jet's path and greatly decelerate the jet head 
and result in a propagation delay which will generally not 
be the same for both jets. We have analyzed the observed 
asymmetry by means of Monte Carlo simulations based on a 
ID propagation model for the jet through the dense medium 
and find that the properties of the ISM can be translated to a 
probability distribution of jet arm length asymmetries. Disc 
density and height are found to have the largest effect on 
the asymmetry for realistic parameter ranges. The shape and 
slope of the Fourier energy spectrum of the ISM density only 
has a small effect on the asymmetry, as has the width of the 
density distribution (related to the turbulent Mach number) 
due to the changes in the filling factor. The hydrodynamic 
simulations showed that asymmetries in general may have 
a wider distribution than expected from the ID model due 
to the complex interaction of the jet and its bow shock with 
gaseous clumps, going much beyond simple energy disposal. 
According to our model, observed asymmetries of medium- 
sized local radio galaxies can be explained by gas masses in 
massive elliptical galaxies of order 10 9 to 10 10 M©. Further- 
more, the hydrodynamic sim ulations provide a theor etical 
basis for the observations of iMcCarthv et al.l (|l99lh that 
emission line nebulae in high redshift radio galaxies are gen- 
erally found to be brighter on the side of the shorter lobe 
due to stronger dissipation and a resulting delay in jet prop- 
agation. This interaction of jets particularly with the cold 
gas phase shows that jet activity may have considerable im- 
pact on the evolution of high redshift galaxies, where star 
formation may be triggered in compressed dense clumps of 
gas, or quenched due to destruction or removal of the cold 
clumps. 



also at larger distances (e.g. stopping cooling flows), radio 
galaxies at high redshift can be expected to show consid- 
erably more interaction with the gas-rich and clumpy en- 
vironment and the cold gas phase. In fact, the clumpy ap - 
pearance of high redshift galaxies (lElmegreen et alJ l2007h 
and clumps for med by Toomre unsta ble discs in cosmologi- 
cal simulations (|Ceverino et al J 1 20 id ) support this scenario. 
And while strong radio source asymmetry is only the mor- 
phological consequence, the actual feedback caused by the 
interaction may be vigorous. 
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6 SUMMARY 

We have performed 3D hydrodynamic simulations of power- 
ful extragalactic jets interacting with a disc-shaped clumpy 
interstellar medium of log-normal density distribution. As 
a first step of analysis, in this paper we have focussed on 
the jet asymmetry and its relation to the inhomogeneous 
ISM. The interaction between the jet and its environment 
does not only occur by direct impact of the jet beam on the 
clouds, but also by a early blast wave phase which consider- 
ably changes the ISM and thereby changes the environment 
through which the beam propagates then later-on. We find 
that the inner region is mostly cleared of dense clumps by 
the jet activity except for some highly compressed structures 
in the disc plane, while the gas at larger radii is only effected 
by the bow shock and the jet backflow impacting on the disc 
vertically. Asymmetries are caused by dense clouds that are 
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